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ABSTRACT 

We examine characteristics of circnmbinary orbits in the context of cnrrent 
planet formation scenarios. Analytical perturbation theory predicts the existence 
of nested circnmbinary orbits that are generalizations of circular paths around a 
single star. These orbits have forced eccentric motion aligned with the binary as 
well as higher frequency oscillations, yet they do not cross, even in the presence of 
massive disks and perturbations from large planets. For this reason, dissipative 
gas and planetesimals can settle onto these “most circular” orbits, facilitating the 
growth of protoplanets. Outside a region close to the binary where orbits are gen¬ 
erally unstable, circnmbinary planets form in much the same way as their cousins 
around a single star. Here, we review the theory and conhrm its predictions with 
a suite of representative simulations. We then consider the circnmbinary planets 
discovered with NASA’s Kepler satellite. These Neptune- and Jupiter-size plan¬ 
ets, or their planetesimal precursors, may have migrated inward to reach their 
observed orbits, since their current positions are outside of unstable zones caused 
by overlapping resonances. In situ formation without migration seems less likely, 
only because the surface density of the protoplanetary disks must be implausibly 
high. Otherwise, the circnmbinary environment is friendly to planet formation, 
and we expect that many Earth-like “Tatooines” will join the growing census of 
circnmbinary planets. 

Subject headings: planetary systems - planets and satellites: formation - planets 
and satellites: dynamical evolution and stability - planet disk interactions - 
binaries: close - stars: individual (Kepler-16) 
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Introduction 


Planet formation is ro bust. Transit detections by t he Kepler satellite (e.g. , 


2011 


2010 


Howard et al.l 120121). radial velocity campa igns fICumming et all 12008 


Zechmeister et al.l 120131: iMavor et al.l 120141 1 , dir ect imaging survey s flMacintosh et al. 


Borucki et al. 


Howard et al. 


2014 lTamurall2014il . and gr avitational len s ing studies flGould et al.ll2010ll suggest that many 


if not all stars host planets (lYoudinl 1201 ll: iDong fc Z 


an 


conhrmed planets, and several thousand candidates (IHan et ahll2014il . 


2013h. To da te, there are over 1500 


The vast majority of the known planets orbit a sin gle star (e.g., iMavor et al.l 12011 


Cassan et al.ll2012l: iBurke et al.ll2014J: iMullallv et ahll2015l) . Their formation in this setting is 


straightforward to describe, even if certain key details are not well und erstood fe.g.. lSafronov 


19691: IWetherilll Il980l: iGoldreich et al.l l2004l: lYoudin fc KenvonI 1201311 . The overall process 
takes many small solid particles of dust and concentrates this mass into a few large objects. 
Coagulation — growth through sticking or merging of planetesimals — is driven by collisions. 
Low relative velocities favor growth, but can also slow it down. Higher velocities can speed 
up the growth rate, but can also lead to destructive collisions. Scattering and gravitational 
interactions pump up relative velocities, while collisional damping and dynamical friction 
slow things down. A balance between these processes enables planets to emerge from the 
dust. 

In additi on to the dozens of pla nets known to orbit one member of a stellar binary (e.g., 
a Cen B; see ISchneider et al.N201lh . a handful of planets are known to orbit both binary 
partners. Kepler-16b - the hrst circumbinary planet disco vered by the Kepl er mission - is a 
Saturn-mass planet at an orbital distance of about 0.7 AU fiDovle et al.ll201l[) . Its binary host 
consists of a 0.7 Mq K-type star and a 0.2 Mq red dwarf at an orbital separation of 0.22 AU. 
Since then, six more “Tatooines” have been reported in the Kepler data set (see Table [H 
below). They are typically Neptune- or Jupiter-size, and all orbit the ir hosts at distance s 
within roughly 1 AU. One binary, Kepler 47, hosts two such planets (lOrosz et al.ll2012al) . 
A few more massive circumbinary planets are known or suspected, but t hese objects orbit 


at mu ch greater distances compared with the binary semimajor axis (e.g., iBeuermann et ah 

201 ih . 


The circumbinary planets that orbit close to their hosts provide unique challenges for 
planet formation theory. The central stars strongly perturb the region around thern , clear - 
ing out orbits to distances of 2-5 times the binary separation flHolman fc WiegertNl999h . 
Similarly, circumstellar disks get eroded from the outside by the binary partner. In either 
case, secu lar excitations from the bina r y potential drive orbital crossings and de s tructive 


collis i ons flMoriwaki fc Nakagawal 120041: iMeschiaril 120121: iPaardekooper et al.l l2012l: iRahkov 


20131: iLines et al.l 1201411 . Thus, planets may not be able to grow near their binary host. 
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To resolve this problem, planets may form at more remote distances, where the time- 
varying part of the binary potential is weak. To arrive at their observed locations, they 
must t hen migrate through the circumbinary gaseous disk or scatter with neighboring gas 
giants dPierens &: Nelson 2008a,b). While the simulations of this process are compelling (e.g., 
Kiev fc Haghighipouiil2014l l. there are uncertainties about starting conditions, typically a set 
of planetary cores placed into a steady state disk. It is unclear if the simultaneous growth 
of planetesimals and dissipation of the disk can conspire to produce cores poised to migrate 
into their observed orbital positions. 

Toward understanding how circumbinary planets form, we re-examine a fundamental 
is sue: the n a ture of pla netesimal orb it s aro und binary stars. Following the approach of 


Lee &: Peald (120061) and iLeung fc Led (120131) . we describe a family of nested, stable cir¬ 
cumbinary orbits that have minimal radial excursions and never intersect. While they are 
not exactly circular, these orbits play the same role as circular paths around a single star. 
Gas and particles can damp to these orbits as they dynamically cool, avoiding the destruc¬ 
tive secular excitations reported in previous work. Thus planetesimals may grow in situ to 
full-fledged planets. 


We organize this paper to provide an introduction to the Lee-Peale-Leung analytical 
theory of circumbinary orbits (§2), followed by numerical examples (§3). We then discuss the 
role these orbits play in planet formation (§4), and close with a comparison to observations 
(§5), along with a summary and predictions of the ideas presented here (§6). 


2. The circumbinary environment 


Planet formation relies on the ability of solid particles — dust, planetesimals, protoplan¬ 
ets — to interact gently. Around a single star, the family of nested (concentric) circular or¬ 
bits offers this possibility. Particles on coplanar circular orbits coexist without any collisions. 
Gravitational interactions among particles (and with coexisting gas) induce random motions 
about these circular orbits, which enables particles to merge into larger objects. Unchecked, 
gravitatio nal interactions grow indehn itely and lead to destructive collisions among parti¬ 
cles fe.g.. iGoldreich fc Tremainel Il978h . Dynamic al cooling (through cqllisional damping, 
dyna mical friction, or gas drag) is essential (e.g., IGoldreich et al.ll2004 lYoudin fc Kenvon 
20131 1. When orbiting particles cool, they need some common set of trajectories on which 
to settle. The family of circular orbits provides this non-intersecting, collisionless haven for 
particles in cold circumstellar disks. 


A central binary dramatically alters these orbital dynamics (e.g., iHolman &: Wiegert 
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19991: iMusielak et al.l 120051: IPichardo et al.l l2005l: IPoolin fc Blundelll I 2 OIII) . For a primary 
and secondary with comparable masses and binary eccentricity, Cbm, satellite orbits are un- 
stable inside of a c r itical radius, Ocrit, which is at least twice the binary separation, Obin- 
Holman fc WiegertI fjl999l) derive an approximation for Ocrit from direct simulation of cir- 
cumbinary particles: 


Ocrit ~ 1.60 + 5.10 Cbin - 2.22 Cbin + 4.12 

M? 




-5.09- 


+ 4.61ebin- 


Mp + Adg 


4.27 Cbin 




Mn + Mg 


( 1 ) 

( 2 ) 


{Mp + MsY + 

where Mp and Mg a re th e masses of the primary and secondary, respectively (see also 
Pichardo et ahl I2005L 120081) . Inside this orbital distance, particles are cleared, creating a 
cavity around the binary. 

Particles beyond the critical distance can be on stable, non-Keplerian orbits. In addition 
to their response to the central mass of the binary, these satellites also experience forced 
motion, driven by the binary’s time-varying potential. This perturbation prevents particles 
from maintaining circular or eccentric orbits. Instead, particles may achieve a “most circular 
orbit,” defined as having the smallest radial excursion about some guiding ce nter, orbiting 
at so r ne constant radius 7?^ and angular speed fig in the plane of the binary flLee fc Peak 
20061: lYondin et ahl 120121) . More generally, eccentric circumbinary orbits may be composed 
of epicyclic motion about Rg, as in the Keplerian case, superimposed on the most circular 
orbit. 


In the rest of this section, we investigate existing analytic theory for circumbinary orbits 
and applications for planet formation. We also include a brief discussion of instabilities and 
resonances, as a prelude to numerical simulations in §3. 


see 


2.1. Analytical theory of circumbinary orbits 

a cen tral binary, we follow the analytic theories of 


To describe satellite orbits abou 


Lee fc Pealel (120061) and iLeung fc Leel (120131), based on the restricted three-body problem 


Szebehelvlll967l: Murray fc Dermottlll999l) . In this framework, a satellite’s position and 


momentum come from equations of motion in the potential of a stellar bina ry. This strat 


1978; 

Marzari & Scholl 

2000; 

Moriwaki & Nakaeawa 

2004; 

Rafikov 

2013) 


e.g., iHeppenheimer 


a satellite’s osculating orbital elements, dehned with resp ect to a central point mass , evolve 
according to an orbit-averaged disturbing function (e.g., Murray fc DermottNl999l) . How¬ 
ever, these elements do not accurately track the orbits of the Kepler circumbinary planets; 
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the binary induces significant motion on dynamical time scales. Although modificatio ns to 
the secular theory can accommodate this extra motion flGeorgakarakos fc Egglll2015h . any 
additional non-gravitational processes like aerodynamic drag or planetesimal collisions are 
more easily described in terms of p ositions and moment a. For these reasons we adopt the 
approach of lLee fc Peald (120061) and iLeung fc Led (120131) . 


The starting point of the Lee-Peale-Leung analysis is the gravitational potential of the 
binary: 


$ = 


GM„ 


GM,. 


+ z^ + R^p + 2RRp cos A0 ^R^ + + R"^, - 2RRs cos Acj) ’ 


( 3 ) 


where G is the gravitational constant and A0 is angle between the secondary and the satellite 
in a reference frame with the binary’s center of mass at the origin. In this frame, the massless 
satellite is at radial position R in the plane of the binary and has altitude 2 ; above this plane. 
The terms Rp and Rg denote the orbital distances of the primary and secondary. 

To make headway, this potential is expanded in terms of the angle cosines , converted 
from p owers (cos^(A0)) to multiple-angle form (cos(fcA0)). For eccentric binaries. iLeung fc Lee 
(I 2 OI 3 I) also expand the potential to first order in the binary eccentricity, eun, using the 
epicyclic approximation to describe the variation in binary separation and phase. Then they 
seek solutions for the excursion of the satellite from a guiding center on a circular orbit of 
radius Rg. The excursions in radial, azimuthal and altitude coordinates are 6 R, 6(j), and 
6 z (which is identically z, since all vertical motions are excursions from the guiding center 
orbiting in plane of the binary). Their solution can be estimated by writing the equations 
of motion and keeping only terms linear in the perturbation coordinates and in the binary 
eccentricity. Between the expansion of the potential and this linearization, the problem 
reduces to the form of a simple, driven harmonic oscillator. 

To follow this prescription, we focus on motion in the binary’s orbital plane. The 
potential is 


<F cos{kA(j)) + Cbin [^Ok cos(A:A0) cos(Gbint) 

-|- 2/c<Fofc sin(/cA0) sin(r2bint)]} [^ = 0] 


fc =0 


( 4 ) 


where hlbin is the mean motion of the binary (Gbin^ = G{Mp -|- Ms)/aun^), and our choice 
of time t fixes the orbital phases of the binary and satellite. The potentials <Fofc are Fourier 
coefficients derived from the expansion of the potential in terms of cos(A0) evaluated at 
Cbin = 0, while the <Fg;i, are those same coefficients giving the first-order terms of a Taylor 
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series in Cbm- Examples of these coefficients are 


"hoo — 


and 


•hoi — — 


GM Gfi 


3{Mp 


Rg: 


■g L 


1 abin^ ^ 9 (M2 + m 2) abin" , 25 (M^ + M^) abin® 


4 Rg 


Gji 

Rp 


Ms) Obin^ 


64 M2 


15 (M" 


Rg 


+ 


256 M4 


Rg 


+ ... 


M 


Rg 


+ 


Mt) 


64 


M4 


Rg 


+ ... 


4*02 — — 


4*03 — — 


$04 — — 


$05 = 


Gj^ 

Rg 

Gfi 

Rg 

Gj^ 

Rg 

Gfx 

Ra 


3 abin" , 5 (M3 + m3) aun^ , 105 (M^ + Mf) 


4R/ 


+ 


16 


5 (Mp - Ms) Obin 


M3 

3 


M 


-Rg 


+ 


R/ 

35 (M 


+ 


512 


Ms) Obin 


M5 

5 


Rg 


+ 


128 


M4 


R/ 


+ ... 


35 (M3 + m 3) abin" , 63 (M^ + Mf) abin'^ 


64 M3 


Rg 


+ 


256 


M5 


Rg 


+ .. 


63 (M^ - M^) abin® 


128 


M4 


Rg 


+ ... 


_ 

^00 ~ " 


Gfi 

Rg 


^01 = - 


^02 


— — 
^03 “ 


g 

G/i 

Rg 

Rg 

G/x 


+ 


16 

Ms) Obin 


+ 


labin^ , 9 (M3 + M3)abin" , 75 (M^ + M^) atin® 
2 

9(Mp 


M3 Rg4 128 M5 
3 75 (M,4 - M^) abin' 


R/ 


+ 


M 


Rg 


+ 


64 


M4 


Rg 


+ ... 


g -‘"'g 

3abin' , 5 (M3 + m 3) abin" , 315 (M^ + M^) abm'^ 


g L 


2 R/ 

15 (Mp 


+ 4' 


M3 


Rg 


+ 


256 


M5 


Rg 


M,)abin" , 175 (Mp" - M,") abin'' 


8 


M 


Rg 


+ 


128 


M4 


Rg 


+ 


^04 


^05 — 


G^ 

Rg 

G/i 


35 (M3 + m 3) abin^ , 189 (M3 + Mf) abin® 


16 


g L 


M3 

r4 


Rg 


+ 


315 (M" - M,") abin' 


256 


M4 


Rg 


128 

+ ... 


M3 




+ 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 

( 10 ) 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


where M = Mp + M^ is the total mass and /i = MpMg/{Mp + Mg) is the reduced mas^. The 
subscripts jk designate that each term is measured in the plane of the binary (j = 0), and 
is the k}'^ harmonic as in Equation (II]). The missing terms are of order (obm/Rg)^, which 
can be as large as a percent in an idealized system, and a fraction of a percent in observed 
binaries. 


^In some previous studies, “/x” is defined as the ra t io of the secondary’s mass to the total mass. Thus, 
our fJ./M is equal to “/x(l — /x)” in Holman fc Wieeert ( IQQOll . for example. 
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From the time-averaged potential ($oo in Equation ([5])), we obtain the angular speed of 
the guiding center, fig. It follows from 




1 (7$, 


00 


Rg dR 


Rg- 


GM 


1 + 


M 


3abin' , 45 (M3 + M3) abin' 


4 R, 


+ 


64 M3 


Rp 


(17) 


The square root of flg^ is the mean motion of the satellite. The satellite’s epicyclic and 
vertical frequencies are 




2 _ 


= Ra 


dQp 


dR 


4Qp^ = 


GM 


Rg 


ly? = 


ld<F 
z dz 


z=0,Rg 


GM 


Rp 


1 + 


1 - 


ii 

M 


3abM , 135 (M3 + m 3) abin" 


4i?g2 


+ 


64 


M3 


Rp 


+ 


ii 

M 


9abM , 225 (M3 + m 3) abin" 


4 i?g2 


+ 


64 


M3 


Rg^ 


+ 


(18) 

(19) 


corresponding to the eccentricity and any motion out of the plane of the binary (for details 
regarding the motion out of the orbital plane, see iLee fc Peald 120061 ). In the limit that the 
binary separation goes to zero, or when the binary mass ratio is extreme, both and z/, 
become the Keplerian mean motion. 


Forced oscillations experienced by a satellite depend on the synodic frequency. 


^syn ^ ^bin 

In general, Wsyn is just the average angular speed of the satellite in a reference frame that 
rotates with the mean motion of the binary. For a circular binary (cbm = 0), the time varying 
force felt by the satellite depends only on this frequency and its harmonics. In the case of 
an eccentric binary, the orbital frequency of the binary also enters into the potential. 



2.1.1. Equations of motion 


Derivatives of the pote ntial yield e q uation s of motion in the excursion coordinates 
{6R,6(f),z). The strategy of Ihee fc Pealel (120061) is to cast these equations in the form of 
a forced harmonic oscillator with natural frequencies Ke (for 6 R and Sf), Ui (for the z coor¬ 
dinate) and driving frequencies involving Dbin and Usyn- Th e solutions in term s of the full 


cylindrical coordinates are (Equations (27), (31) and (35) in lLeung fc Leell2013l) : 


7^(^) Rg \ 1 6free OOsi^K^t ^ ^ G}^ COs(A)Cc)gynt) (^^) 

I k=l 

~ oo 

Cl COs(Dbin^) + ^ Cf COs{kUsynt Dbin^) + Cjl COs{kUsynt - ^hint) 


^bin 


k=l 

































oo 


Du 


(j){t) = f2g t Sm{Ket + V’e) + / T- 

Kg kUgyyi 


sin(fcWsynf) 


( 22 ) 


“l“^bi] 


■^0 


^hi 


COS 


bin 


/O sin(fcWsy„t + l^binb , sin(A;a;synt - ^binb 

(.“bintj + ^ -;-—-H 


k=l 


kuJsyn + D, bin 


ku, 


'syn ^bin 


z{t) = iRg cos(i/it + 'Ipi), 


(23) 


where Cfree is the “free” eccentricity, i is the inclination, and the phase angles i/'e and xjjt are 
constants. The coefficients are 


C, = 


f^e — 
Uq — 


1 


R^{kI - 

1 

Rginl - ribin^ 




Ofc 




g'nofc 


diti kigCO^y]^ 




Rs 


dR 


Re 


ci 

fkK 

Dk 

= 2Ck + 


= 


= + 


±k 




Ofe 


dR 2 dR 


Id^ _ fcl]g(±2fc$ofe - ^^) 
-R(fccJsyn i ^bin) 


Re 


$ 


Ofc 


2 

^2gCt7gyj^ 


Rr 


k{±2k^ok - 
_2Rg r2g(fccJsyn i fibin) 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


Rs 


(Equations (28-30) and (32-34) in Leung fc Lee 2013 )1^ 


2.1.2. Components of the orbital motion 


From the lT;enng fc T^eel (120131 1 solutions, we see that circumbinary orbits may be broken 
up into independent modes. The hrst mode is forced motion, prescribed by the characteristics 
of the binary and the orbital distance of the satellite from the center of mass. The second 
mode is “free” motion, fully analogous to eccentricity and inclination of orbits around a 
single central mass, except that the epicyclic excursions are relative to an orbiting reference 
frame locked into forced motion, as opposed to a circular guiding center. 


^The coefficients C and D are equal to the corresponding C and D in lLeung fc Leel (120131) divided by the 
binary eccentricity. For example, C'^!’ = Cf /cbin and Cg = Co/ebin- Our choice allows the dependence on 
binary eccentricity to appear explicitly in the solutions of the excursion variables. 
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The forced motion itself can be broken down into parts, including fast, driven oscilla¬ 
tions at the synodic frequency and the binary’s orbital frequency, plus the slower epicyclic 
oscillations operating at the orbital frequency, fig. The higher frequency contributions have 
radial excursion amplitudes that scale as 


^2, Ci^ !> X R. 


A 

M 


®bin 

Rp 


Rp 


(30) 


to leading order in a^nlRg le.g.. Iheung fc heel l20 1311 . The slower epicyclic motion is asso¬ 


ciated with the forced eccentric orbit (iHeppenheimerlllQTSl ). Its contribution to the radial 
excursions is typically much larger, scaling as 


^force^^g Rg 


5 - M, 


M 


"^bin^bin* 


(31) 


When the binary eccentricity is large, we rewrite the solutions to highlight the role of 
the forced eccentricity. For example, the radial coordinate becomes 


R(t^ Rg [l Cfree COsi^K^t "fpe) ^force COs(flgt) “1“-] 


(32) 


This form of the solution is familiar from secular perturbation theory (IMurrav fc Dermott 


19991). 


An important feature of the forced eccentric orbit is that it does not process. Its argu¬ 
ment of periastron is locked in line with the argument of periastron for the binary. All other 
components of the forced motion are also synchronized to the binary’s orbital frequency, the 
synodic frequency or their harmonics. In this way, the orbital paths of satellites experiencing 
only forced motion make a family of nested orbits that never intersect. 


In contrast to the forced eccentricity, the free eccentricity and inclination are both 
associated with precession. This effect is a direct result of the fact that the time-averaged 
potential around a binary does not fall off as l/i? — t he binary’s mass averaged over its orbit 
is ak in to an oblate spheroid flHeppenheimeii Il978l : iMurrav fc DermottI Il999l: iLee fc Peak 
2006h . The precession rates of the periastron and ascending node are 


VJ = 


Q 


node — 


3 ttbin /X 

'CM/Rg^ 

(33) 

® A R,'^ M 


(34) 


and thus are approximately equal and opposite. 
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2.2. Incorporating a gaseons disk 


The orbit solutions in Equations fl?I]) - fl2^ apply directly to non-interacting particles in 
a disk with negligible mass. Now we consider how gravity, pressure support, and aerodynamic 
drag from the gas disk modihes these orbits. To make headway, we assume that outside of 
the critical radius where orbits are unstable, the binary’s gravitational potential changes the 
fluid flow in a disk from circular orbits to most circular paths, similar to the satellite orbits 
described above, but dependent on the physics of the disk. Our premise is that fluid flow 


of free eccentric orbits (see discussions in 

Lin & Pringle 

.976. Sver & Clarke 

1992, 

Ogilvie 

2001 , 

Ogilvie & Barker 

2014 

, and 

Barker & Ogilvie 

201 ^ 

1). Thus nested, non-intersecting. 


most circular streamlines are a natural extension to circular flows around a point mass. 

Continuing with this picture, we make the simplifying assumption that to a good ap¬ 
proximation the disk potential and gas pressure are axisymmetric. If gas fluid elements 
travel on most circular orbits, then radial and azimuthal variations in disk properties must 
exist. However, they may be small; in the Kepler circumbinary systems, the forced ec¬ 
centricities are Cforce ~ 0.002-0.044 (see Analytical models conhrm that close to the 
binary, the gravitational effects stemming from disk eccentricity are small compared with 
the stars’ influence on circumbinary orbits. In this way, we build on previous descriptions of 
circumbinary disks that assert strict axisy mmetry in the disk by having fluid element s travel 
in pressure-supported circular orbits fe.g.. iMarzari fc Schollll2000l : IScholl et al.ll2007t) . 


In the limit of a tenuous gas, this simple picture yields reasonable results. As pressure 
and disk gravity diminish, gas molecules damp to the same set of most circular orbits derived 
for satellites in ^2.11 However, the model does not incorporate hydrodynamical effects, 
including viscosity, turbulence and gravitational instab ilities, that may signihcantly modify 


disk structure fe.g.. iPelupessv fc Portegies Zwartll2013|) . With this acknowledgement of the 


limitations of model, we proceed to consider circumbinary orbits. 


2.2.1. Orbits in the presence of a massive disk 

The gravity of a massive circumbinary disk modihes the orbits of satellites around a 
binary star. An axisymmetric disk, with a gravitational potential ^d{R) in the orbital plane 
of the binary, modihes the mean motion and the epicyclic frequencies of a circumbinary 
satellite. These properties of the satellite are related to derivatives of the time-average 
potential; we derive them by making the substitution 


d>oo ^>00 + 


(35) 
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in Equations ffTTI) and ffTSl) . These changes to fig and Ke cause only minor adjustments to 
most terms in the orbit solutions except for the term associated with forced eccentricity: 


cr = 


{kI - fig^) L 


d 

dR 


$ 


01 


+ 

2 


01 


+ ^ (at’oi - 4ra) 


(36) 


Re 


With fig and in the denominator, the magnitude of the forced eccentricity now has a wide 
range of va lues tha t depe nd on the form of the disk potential. In plausible astrophysical 
conditions, iRafikovI (120131 ) demonstrates that the denominator can be negative and large 
compared to a disk-free system. The magnitude of the forced eccentricity then becomes 
small; the forced eccentric orbit is anti-aligned with the binary. 


When the denomina tor in Equation fl36|) goes to zero, satellite orbits experience a secular 
resonance fjR,afikovll2ni3l ). Physically, this condition occurs when orbits with free eccentricity 
do not process; contributions to the apsidal precession from the disk and from the binary 
are equal and opposite. To quantify this behavior, we choose a disk surface density of 


E,(a) = 

(Jj 


(37) 


with Sg o = 2000 g/cw? and Oq = 1 AU, typical of observed gas disks le.g.. lDent et al.ll2013l) . 
The disk potential is 


= 27iGT,g^o aolog(a/ao) 


(38) 


flBromlev &: KenvonI l2011ai Appendix A). The corresponding apsidal precession rate from 
the disk is 


'^d 


d 






,d^d 

dR 


TT 


Re 


r2gi?g 


GS(i?g). 


(39) 


Setting the magnitude of this expression equal to the apsidal precession rate from the binary 
(Equation ([33])) yields the radial position of the resonance. 

For the Kepler circumbinary planets, the resonance lies at an orbital dista nce of roughly 
2 AU. This distance is bet ween the orbits of the know n planets (lR,afikovll20131 ) and the likely 
position of the snow line flKennedv fc Kenvonll2008h . As the disk dissipates, however, the 
resonance sweeps outwards past the snow line. This motion may have interesting implications 
for the formation and inward migration of gas giantsj^ 


^It is beyond our scope to treat the orbital dynamics at the resonance in detail. Aside from its impact 
on the orbits of solids, it is unclear how the gas reacts. Within the Lee-Peale-Leung theory, the evolution 
equation for the Ci mode on resonance is an undamped, driven harmonic oscillator with a natural frequency 
of ftg. This solution predicts a linear growth time scale for ec centricity of 10 ^-10^ years for the Kepler 
circumbinary planetary systems, which may be an overestimate (lMeschiarill20l4 . Figure 6). 
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An extension to the axisymmetric form for the disk potential would account for the 
non-axisymmetric ebb and flow of fluid elements on most circular orbits. A starting point 
is to approxi mate the streamlines with the forced eccentric orbits from secular perturbation 
theory (e.g., ISilsbee fc Rafikovll2015blJal) . The gravitational held of the eccentric disk then 


affects the forced eccent ricity, but not the apsidal alignment (see Equations (19)-(21) in 


Silsbee fc RahkovI l2015al) . We could then calculate a self-consistent value for Cforce where 


the orbits of the huid elements follow paths that they help to generate. Further extensions 
to the theory of eccentric orbits would include the effects of geometric compression of huid 
density from the co ntinuity equation and related hydrodynarnical e hects for how along most 
circular paths (e.g., lOgilvie fc Barkerll2014J: [Barker fc Qgilviel 120141) . 


S.S.S. Orbits in the presence of a pressurized disk 

The pressure in protoplanetary gas disks modihes the orbit of any huid element. In a 
simple axisymmetric circumstellar disk where t he radial pressure gradient is positive, pres¬ 


sure support leads to su b-Keplerian orbits (e.g.. IWeidenschilling||l977al: iBirnstiel et al.ll2010 


Chiang fc Youdinll2010[) . The mean orbital speed of the gas relative to a circular Keplerian 
orbit is 

|An| r/i?gf2g, (40) 

where rj ~ 10“^ fe.g.. IWeidenschilling||l977al) . In a circumbinary disk, we estimate the ehect 
of pressure by assuming that the pressure gradient is radial, and treating 77 as a constant, 
independent of orbital position, at least in some local region of the disk. We then may let 
pressure support appear in the equations of motion as a change in the total mass of the 
binary: 

*^*00 (1 — 2?7)<hoo, (41) 

leading to small modihcations in quantities including fig and Thus, gas pressure increases 
the forced eccentricity by a factor of 1 -|- 277 compared to an unpressurized disk. 

In this approximation, t h e gas and solids not coupled to it are on d istinct most circular 
orbits f Thebault et ah 2006 : Meschiari 2014 : Silsbee fc Rahkov 2015a '). Nonetheless, the 


forced eccentric orbits of disks with and without pressure support are apsidally aligned. 
Furthermore the diherences between forced eccentric paths of gas and the uncoupled solids 
are small. They lead to relative radial speeds between gas and these solids of 


Auj. ~ 


(42) 


which, for anticipated values of Cforce for the Kepler circumbinary planets, is roughly two 
orders of magnitude smaller than the “headwind” felt by the solids as they plow through 
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the gas. In other words, the most circular paths associated with pressure-supported gas 
and uncoupled solids are sufficiently similar that relative velocities arise predominantly from 
their azimuthal motion, just as in circumstellar disks. 


2.2.3. Gas drag 


When the gaseous disk has internal pressure, material orbits the central object more 


slowly than the solids . The solids then experience aerodynamic drag flAdachi et al.Nl976 


Weidenschi]]ingill977aj) . Particles smaller than about a centimeter are fully entrained in the 
gas. Planetesimals with radii of 1 km or more barely feel the gas. Both sets of particles 
follow distinct most circular orbits. For intermediate particle sizes, aerodynamic drag slows 
azimuthal speeds. Without the benefit of radial pressure support, these objects spiral inward 


(Adachi et al. 

1976; 

Weidenschilline 

1977a; Weidenschilling & Cuzzi 

1993: Chiane & Youdin 

2010 ; 

Youdin & Kenvon 

2013 

)• 


To describe the dynamics of a particle experiencing gas drag, we take advantage of the 
small difference between the forced eccentricities of a pressurized gas disk and satellites that 
orbit in the absence of gas drag. Approximating the effect of gas as a constant azimuthal 
headwind, the equations of motion resolve to a most circ ular orbit with some inte rmediate 
forced eccentricity and a small amount of radial drift fe.g.. lYoudin fc Kenvonll2013h . Impor¬ 
tantly, if the epicyclic motion of the entrained particles and large solids are apsidally aligned, 
then so are the epicyclic orbits of these intermediate-size bodies. We confirm below using 
numerical tests that this description is reasonable (^. 


2.3. Significance for circnmbinary planet formation 

Here, we highlight several features of the orbit solutions described in 1 12.11 and modifi¬ 
cations arising from the presence of a gas disk fl l2.2p : 

• The forced epicyclic motion (the (Pf term) is synchronized with the guiding center. 
The addition of an axisymmetric potential, whether it describes the effects of a massive 
disk or mimics the behavior of radial pressure, preserves this relationship. In general, 
forced eccentric orbits remain apsidally aligned with the binary. In a disk free of gas 
pressure, particles on these paths never collide. 

• Pressure support changes the most circular paths of the gas streamlines relative to the 
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orbits of solid particles not susceptible to aerodynamic drag. This behavior generates 
relative epicyclic motion between the gas and the solids. If the forced eccentric orbits 
are apsidally aligned with the binary, however, then the relative speeds induced by the 
differences in the most circular paths are much smaller that the usual headwind felt 
by the solids plowing through the gas. 


The forced epicyclic motion is coupled in its alignment with the binary. Thus if the 
binary itself precesses slowly, then the satellite’s argument of periapse also processes. To 
visualize this point, we note that a slowly processing binary orbit has slig htly different 
radial and azimuthal frequencies. To the solutions of lTjOnug fc T^eel (120 131) , we thus add 
the precession rate Whin to the binary’s orbital frequency, hlbin, without modifying the 
synodic frequency (Ugyn between satellite and binary. The terms most greatly affected 
by this small change are the ones associated with the forced eccentricity (e.g., with 
Cl). The time dependence in these terms then transforms as 


cos(a;synt - flhint) -t cos((f2g + d7bin)t). 


(43) 


Thus the satellite’s forced epicyclic motion precesses with the binary, remaining apsi¬ 
dally aligned (in secular resonance). 


These properties may be essential for circumbinary planet formation. They suggest that 
most circular orbits remain nested and non-intersecting even in the presence of a massive 
disk or if the binary precesses due to interactions with a massive disk or a distant plane¬ 
tary/stellar perturber. Gas dynamics add complications, although if a protoplanetary gas 
disk is nearly axisymmetric and streamlines follow most circular orbits, then the presence 
of the binary induces only small additional relative speeds compared to the circumbinary 


lar planet formation (e.g.. 

Youdin & Chiane 

200 

4; 

Birnstiel et al.l 

2010 ; 

Ghiane & Youdin 

2010 ; 

Windmark et ah 

2012 

; Garaud et ah 

20131 

Youdin & Kenvon 

2013 

carry over to the 


circumbinary environment. 

Other studies of circumbinary planet form ation include disks that are exactly axisym- 


metric and circular in their orbital flow (e.g., IScholl et al.l l2007l) or that are neither ax 


isymmetric or apsidally aligned with the binary flSilsbee Sz Rafikovll2015a|) . In these cases. 


aerodynamic drag on particles leads them to achieve orbits that depend on their physi¬ 
cal size. In misaligned disks, for example, the magnitude of the force d eccentricity and 


2015a 

). In either case, p 

(Marzari & Scholl 

2000l 
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We assume that most circular orbits, aligned with the binary, are a good hrst approximation 
to the paths of fluid elements. In any event, as the gas dissipates, we expect that the fluid 
elements and solids are likely to settle on the same set of most circular orbits. 


2.4. Summary of the analytical theory 


The analytic theory of circumbinary orbits predicts a family of nested most circular 
paths. Satellites orbiting with these paths (i) make well-defined minimal radial excursions 
and (ii) never collide. Thus, these paths dehne “dynamically cold” orbits in exactly the 
same way as circular orbits around single stars. As in the Keplerian case, satellites may have 
additional “free” eccentricity and inclination to describe motion about these paths. For any 
combination of free and forced eccentricity, the dynamics of planetesimals as they stir or 
damp each other takes place in the frame of these most circular paths. 


While the physics of gaseous protoplanetary disks is uncertain, gas fluid elements can fol¬ 
low most circular paths, even when the disk mass and pressure support are signihcant. These 
orbits serve as reference frames for local hydro- and aerodynamics. We expect that gas settles 
to these orbits, since they allow for streamlines to be nested and non-crossing with minimal 
radial excursi ons, despite any forced eccentric motion. Studies of eccentri c gas disks around 
a sin gle star f Sver fc Clarke 1992 : Qgilvie 2001 : Qgilvie &: Barker 2014; Barker fc Qgilvie 


20141 ) highlight potential differences between the hydrodynamics of circular flows and eccen¬ 


tric ones, including ma ss conservation along eccentric streamlines fjOgilvie fc Barker! I2014J: 


Barker fc Qgilviel 120141). In circum stellar disks, gas might find stable orbits with some free 


eccentricity fjSver fc Clarkd Il992 ) . Ho wever, it seems more likely that apsidal precession 
circularizes these orbits (jOgilviel 120011) . For circumbinary disks, we expect that apsidal 
precession leads to fluid flow along most circular paths. 


2.5. Limitations of the theory: resonances and chaos 


Limitations of the theory stem from its perturbative approach, particularly that it is 
linear in eccentricities. Aside from singularities in the coefficients in the solutions (e.g., Ck, 

by a massive dis k flRafikovi l2013l. see 1 12.2. 11 above 


^ ) at the 1:1, 1:2 and 1:3 comm ensurabilities — along with the possible resonance induced 


resonant effects fjWisdom 


1982 


1983 


^2.2.11 above 

— t 

Lecar et ah 

2001 ) 


These phenomena would arise 


in the equations of motion if the satellite’s eccentric motion were explicitly includ ed in the 
calculation of the force. Then, solutions may become chaotic and possibly unstable flWisdom 
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iQsnh . 


The formula for arrit provides a starting point for investigating instability (Equation [T| 
Holman &: Wiegertlll999[) . In this approximation, the location of unstable orbits close to the 
binary varies smoothly with the mass ratio and binary eccentric i ty. However, locating all 
of the unstable or 


Dits is more complicated (IMusielak et ahl l2005l: iDoolin fc Blundell 


Chavez et ah boiS h Overlapping resonance conditions generate instability f Chirikov 


2011 


1959 


19791: IWisdomI Il980l) : thus, unstable orbits can exist in narrow, isolated ranges of orbital 
distance a beyond Ocrit- hi the binary Kepler-16 (Tabled]), the 5:1 resonance, located out- 
side of arriti is unstabl e. Although it is weaker, the 6:1 resonance is also unstable (e.g., 
Pierens fc Nelsonll2007l) . In b etween these resonances, orbits are stable: Kepl er-16b resides 
on one of these stable orbits flPopova fc Shevchenko! l2013l: IChavez et alJl2015l) . Beyond the 
6:1 resonances, all orbits are stable for small planetary eccentricities. When we discuss 
formation mechanisms for circumbinary planets, we return to this issue (]|5|). 


3. Numerical simulations 

The analytic theory for circumbinary orbits is hrst-order accurate in the binary eccen¬ 
tricity and is derived from equations of motion linearized in the excursions away from the 
guiding center. Despite these limitations, it compares well in numerical experiments even 
when the binary eccentricity is moderately large. We summarize these experiments in this 
section. 


3.1. Code description 


We apply our planet formation code Orchestra flBromlev &: Kenvonll2006l. l2011al) to 
compare with the theoretical results in Orchestra is a hybrid u-body-coagulation code 
for tracking the emergence of individual planets (n-bodies) from a sea of smaller particles 
that can be characterized statistically (the coagulation “grid” with bins for radial position 
and particle mass). In a standard hybrid calculation, small particles within a set of concentric 

ensity and initial e and i relative to a circular 


orbit (see also 

Safronov 

19691: 

Lissaue 

1987 

: Suaute et ah 

1991 

: Wetherill & Stewart 

1993; 

Weidenschilline 

1989; 

Kenvon & Luu : 

1998; I 

Cenvon & Bromlev 

20081). Initiallv, there are no 


u-bodies. As massive objects evolve in the grid, the most massive are “promoted” into the 
u-body part of the code. The subsequent evolution of the n-bodies and the grid are linked 
together, enabling the simultaneous tracking of gravitational interactions and collisions that 
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lead to accretion, merging and fragmentation (e.g., iKenvon &: Bromlevll2006L |2010[ l2014l) 


Orchestra has other capabilities which allow us to consider a variety of problems in 
planet formation and evolution. It can track a swarm of massive or massless tracer particles - 
usually sampled from the orbital distribution of the coagulation particles - to mediate inter¬ 
actions that involve resonances and migration flBromlev fc Kenvonll2011br). The code also in- 
cludes interactions with a massive gas disk, both through gravity (Bromlev fc Kenvon 2011a ) 
and aerodynamic drag fjWeidenschillind Il977al: iKenvon fc Bromlevi 1200211 . Tests of these 
and other elem ents of the code, including our choice of time integ rator — either symplectic 


1 rosmaa ryyu i or aa< 

— are summarized in 

iprive JAicna 

Kenvon & E 

rason ex 

Iromlev 

rrapoi 

(2001) 

arion at o -oraer i ai 

. Bromlev & Kenvon 

omiey 

(2006) 

(V: jxenvoni zuud i 

. Bromlev & Kenvon 

(2011al). and Kenvon & Bromlev 

(20151). 


Here we use the n-body component of the code with tracer particles to calculate orbits 
around a binary. The primary and secondary stars are n-bodies, evolved with the 6*^-order 
symplectic integrator. Energy errors are better than one part in 10^°. In some runs (see 
below) we include a third, Jupiter-mass n-body as a perturber. In others we include a massive 
gas disk, choosing a surface density E in Equation (157)) and the corresponding gravitational 
potential in Equation fl38|) to use in the equations of motion for the tracers. This approach 
allows us to test the analytic theory in 1|2] using orbit solutions derived numerically at high 
accuracy. 


3.2. Simulation results 


Simple n-body experiments with Orchestra allow us to test several key features and 
predictions of t he analytic theor y (§2). For these studies, we use the Kepler-16 system 
as an exampl e JPoyle et 2011 ). adopting binary parameters in Table [1] as estimated by 


Leung fc Leel f)2013l) . 


We begin with several illustrations of circumbinary orbits with no other massive per- 
turbers. Our first example focuses on most circular orbits, showing both high-frequency 
oscillations and epicyclic motion. In a second example, we consider the difference between 
the time evolution of particles on most circular and ‘initially’ circular orbits around a ec¬ 
centric binary. Particles on most circular orbits do not process. Particles which start out 
on geometrically circular orbits have some free eccentricity; this component of the orbit pro¬ 
cesses. By selecting an initially circular orbit with equal parts of free and forced eccentricity, 
we show how the precession of the free eccentric orbit modihes trajectories around the binary. 


This second example identihes issues with previous n-body studies of circumbinary 
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plane t formation (e.g., Moriwaki fc Nakagawa 2004 : Scholl et al.lboOT : Meschiari 2012 : Paardekooper et al 


20121) . In these analyses, initially geometrically circnlar orbits process, cansing high relative 
collision velocities that can inhibit the growth of planetesimals into planets. Most-circnlar 
orbits have smaller radial excnrsions and provide a calmer frame of reference for planet 
formation. 


We then explore the impact of a massive disk, gas drag, and a Jnpiter-mass pertnrber. 
Onr goal is to investigate whether most circnlar orbits remain “most circnlar” in the context 
of a protoplanetary disk aronnd a binary star. 


• A sequence of most circular orbits. Fignre[T] shows the resnlt of test particles on most 
circnlar paths near the orbital distance of Kepler-16b. In these orbits, most of the 
motion comes from the forced eccentricity. As snggested by the Fignre, the orbits are 
all nested and do not intersect. 


• A most circular orbit in detail. Fignre [2] contains radial excnrsions of a satellite of 
Kepler-16 where the binary’s eccentricity is rednced by a factor of ten. With this low 
eccentricity, both the high-freqnency driven oscillations and the forced eccentric orbit 
are apparent. This Fignre shows a direct comparison to theory (Eqnation fl^T]) h 


• “Circular” versus Most-circular. Satellites can be lannched on paths that are initially 
more circnlar (smaller radial excnrsions) than the most circnlar orbits discnssed here. 
However, these satellites do not remain on circnlar paths as their orbits evolve over 
time. Fignre [3] illnstrates this point. For a central binary with non-zero e, an orbit with 
eqnal parts of free and forced eccentricity can be initialized on a pnrely circnlar orbit 
abont the binary center-of-mass. At the start (time = 0 in the hgnre), this orbit has 
very little epicyclic motion (e.g., dR is close to zero). Over time, the free eccentric orbit 
processes (Eqnation (133D ): the forced eccentric orbit does not. Hence the two modes 
drift in phase, cansing the beat pattern in the radial excnrsion shown in the Fignre. 
Thns, orbits can be geometrically “circnlar” bnt only temporarily. Most-circnlar orbits 
have the smallest radial excnrsions over the long term. 


Previons stndies, inclnding onr own flKenvon fc Bromlevi 1201411 . describe simulations 
to track particle dynamics around binaries. If particles are set up initially on geo¬ 
metrically circular orbits (where free and forced eccentricities cancel at t = 0), they 

“secular excitations” 


time until e = 2efr,rrp 

(Moriwaki &: 

N'akaeawa 

2004; 

Scholl et al. 

2007; 

Meschiari 

2012 ; 

Paardekooper et al. 

2 OI 2 I: 

Rahkov 

20131. However, the particles are never actuallv 


stirred by the binary. They simply experience the independent free and forced modes 
of epicyclic oscillation, acting in concert but not in phase. 
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• Most-circular paths around a high-eccentricity binary. In Figure IU the eccentricity of 
the Kepler-16 binary is increased to cun = 0.5. The satellite is at 0.9 AU, just outside 
the critical radius for stable orbits for this choice of Cbin. The Figure shows samples of 
the planet’s position in the plane of the binary, taken over the course of 10^ satellite 
orbits (blue points in the Figure). The samples reveal the forced eccentricity and 
demonstrate that the periastron of the satellite is hxed and aligned with the binary. 
This numerical test is signihcant: the analytic theory described here is linear, but the 
n-body experiment is not. The orbital alignment — a prediction of the linear analytic 
theory — holds in the non-linear case. 


Orbits in a massive disk. Figure H] also demonstrates that most circular orbits do not 
process even in the presence of a massive disk. We use a potential as in Equation flHTll 
with a surface density of 2000 g/cm^ at 1 AU. Around a single star with Kepler-16’s 
mass, the disk causes rapid apsidal precession of a satellite at 0.9 AU, do ~ 0.005 yr“^ 
flRahkovI l2013l ). However, around the binary, the satellite’s forced eccentric orbit is 
immune to precession from the disk, as we expect from perturbation theory. The 
example in the Figure demonstrates that this prediction extends to moderately large 
values of the binary’s eccentricity. 


• Gas drag. Aside from changing the gravitational potential, gas modihes the orbits 
of small solids through aerodynamic drag. Figure [5] shows a particle orbit in a sub- 
Keplerian (“pressurized”) gas disk whose mean azimuthal speed is reduced by a factor 
of 1 — 7] relative to a guiding center orbit in the disk’s absence {p = 0.001; Equa¬ 
tion 0401) h For the purpose of this illustration, we do not include the disk’s gravity. 
The drag force is proportional to the particle’s speed in the local most circular refer¬ 
ence frame of the gas. The magnitude of the force is sufficient to cause the particle to 
inspiral within a few hundred orbital periods. The Figure illustrates that even in the 
presence of gas drag, the particle orbits remain aligned with the binary. 


• Effects of stirring: external time-dependent perturbations. Collisional damping drives 
particles toward most circular orbits because particles coexist on these orbits without 
collisions. Conversely, distant massive perturbers (i.e., planets) gravitationally stir 
particles, driving them away from most circular paths. From the analytic theory, we 
expect stirring to behave the same way around a binary as it does around a single star. 
Figure [6] shows simulation data from three scenarios: (i) a satellite like Kepler-16b; (ii) 
a satellite together with a Jupiter mass-body on a circular orbit at 2 AU; and (iii) a 
satellite and a Jupiter orbiting a central point mass. Perturbations of the more distant 
planet affect the motion of satellites in circumstellar and circumbinary cases in much 
the same way, by generating nearly identical free eccentricity. 
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• Binary precession. When an eccentric binary precesses, we expect the satellite’s forced 
eccentric orbit apsidally precesses at the same rate. Figure [7] illustrates this behav¬ 
ior. In a reference frame rotating with a processing eccentric (cbin = 0.5) binary, the 
satellite’s orbit shows no precession due to a Jupiter mass perturber at 2 AU. 


3.3. Summary of the numerical studies 

To summarize, this set of numerical experiments illustrates that the main predictions 
of the linear analytic theory are conhrmed in the non-linear regime. A key result is that 
particles initialized on geometrically circular Keplerian orbits have a free eccentricity, which 
leads to precession and possibly high-velocity collisions among particles on adjacent orbits. 
In contrast, most circular orbits remain nested and never cross, even in the presence of a 
massive disk or a gas giant perturber. Our conclusion is that the standard initial conditions 
of particles on geometrically circular orbits are not realistic in a dynamically cool planetary 
disk. The nested most circular orbits are a better starting point. In the next section, we 
examine these conclusions in the context of planet formation theory. 


4. Circumbinary planet formation 

In the standard theory of star and planet formation, a rotating m olecular cloud of gas 


1981; 

Terebev et ah 

1984; 

Yorke et ah 1993). Disk material is on i 

5 ^ 

learly circular orbits 

(e-g-, 

Weidenschilling 

1977 

a): viscous shear transports mass inward 

and angular momen- 
3 led to and flow with 

turn c 

the gc 

utward 

Lvnden-Bel 

& Pringle 1974). Small solids are well-com 

IS (e.g.. 

A.dachi et ah 

1976: Weidenschilling 1977a; Rahkov 2004: 

Birnstiel et ah 2010: 

Ghiang & Youdin 20101). Larger solids decouple from the gas and follow Keplerian orbits 


about the central star. 

Among the decoupled solids, various dynamical processes tend to circularize their orbits 
around the central star. Smaller particles feel a strong headwind and are dragged towards the 
cent ral star. Although larger particles feel less drag, the gas effic i ently damps their orbits (se e 


also Safronov 1969 


Greenberg et al.l Il978l: ISpaute et al.l Il99ll: IWetherill fc StewartI Il993l) . 


Along with gas dra g, collisional damping and dynamical friction drive particles towards 


circular orbits (e.g.. iHornung et al.l Il985l: IWetherill fc StewartI Il993l: iKenvon &: Luul Il998 


Goldreich et al.ll2004|) . As a result, most particles ex perience small relative co llision velocities 
which encourages growth through mergers (see also lYondin fc Kenvonll2013l) . 
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When the central protostar is a close binary, we expect a similar evolntion. Infalling 
gas forms a circumbinary disk where the gas and dust generally follow most circular paths 
around the binary. Aside from an inner gap in the disk roughly at Ocrit, the structure of 
the disk orbiting a circular binary is fairly sim i lar to a circumstellar disk a round a single 
star fe.g.. lLin fc Papaloizoulll979l: lPringlelll99ll: lArtvmowicz fc Lubowill996h . For eccentric 
binaries, disk material appe ars to follow most circu l ar orbits with forced eccentricity, d riven 
by the central binary fe.g.. iPierens fc NelsonI 120071: iPelupessv fc Portegies Zwartll2013|) . As 
long as the disk is not dominated by non-axisymmetric structure (e.g., spiral density waves in 
a massive disk), we expect pressure and viscosity to induce a smaller inward drift of t he gas 
and small particles relative to most circular orbits (see also iPichardo et ahll2005L l2008l) . The 
gas attempts to circularize the orbits of larger particles onto most circular orbits. Collisional 
damping and dynamical friction also damp the orbits. Thus, large particles end up on most 
circular orbits with a small amount of free eccentricity with a magnitude similar to the 
eccentricity of particles in disks around a single star. 


Achieving these configurations is a natural long-term outcome for dissipative disks. How¬ 
ever, it is worth considering how quickly the gas, dust, and larger solid particles circularize. 
We then compare these time scales with time scales for other processes such as precession 
of the binary. We expect that particles settle onto most circular orbits when the precession 
time scales are long relative to circularization time scales. Here, we consider several basic 
time scales, treating gas-dominated and particle disks as separate cases due to differences in 
their surface densities and in the physical processes that drive them. 


4.1. Gas disks 


We start by con sidering a massive circumbinary disk to model primordial gas in pro¬ 
toplanetary systems fjWeidenschi]]indll977bl: lHavashilll98ll) . As before, we choose a surface 
density 

-1 




Oq 


where q = 2000 g/cm^, and Oq = 1 AU. The disk has a vertical scale height 


h{a) = hg^o 


ao 


(44) 


(45) 


where hgfi/ao ~ 0.02 and q = 9/7 fIChiang fc Goldreichlll997h . 
Damping times in this disk derive from the sound speed. 


Cs(a) ~ h{a)v^/a, 


(46) 
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where is the orbital velocity at distance a. The time scale for vertical structures to reac h 
hydrostatic equilibrium is h/c^, close to the dynamical time flLightmanlll974J: |Pringlelll98lh . 
Pressure damping of features along an orbital path require at least 


Tg,damp ^ o/Cs ~ ^ 7 ^ ~ 


■ a ■ 

17/14 

■ M ' 

LiauJ 


_lMo_ 


1 - 1/2 


yr. 


(47) 


where M is the t otal central rn ass. Radial structures may dissipate more slowly, on a diffusive 
time scale (e.g., IOgilviell200lh 


^ p,damp 


< ^2.5x10=^ 


■ a ■ 

13/14 

■ M ' 

ll AuJ 


_1 Mq_ 


- 1/2 


yr. 


(48) 


We anticipate that disks with free eccentricity damp to most circular orbits on time scales 
in this range. 

By comparison, the time scale for binary precession, a result of the gravitational inter¬ 
action between the gas disk and the binary, is long. If the inn er edge of the disk is at orb ital 
distance Oin, then the binary’s apsidal precession rate is fe.g.. IWardlll98ll: iRafikovI 120131) : 


^bin 


0.5 


7rGSg(flin) ^®bin 

f^bin® in 


(49) 


The precession time is then 


Tbin-pre ~ 2.0 X 10^ 


-1 

M 

O 

-1 

Obin 

1/2 

flin 

2000 g/cm^ 


Lo.2AuJ 


_2.5abin. 


yr [gas disk], (50) 


where we choose to set Oin to be the innermost stable orbit at Ocrit ~ 2.5abin, for the values 
of the parameters in the angular brackets. 

Thus, binary precession is slow compared with dynamical times. Precession is a per¬ 
turbative effect. We therefore expect that orbits in the disk remain apsidally aligned with 
the binary. We conclude that gas disks damp quickly to most circular orbits, so long as the 
physics (gravity, hydrody namics) enables the disk to be axis ymmetric when averaged over 
these orbits (although see iPelupessv fc Portegies Zwartll2013l) . 


4.2. Particle disks 


As the gaseous disk evolves, solids particles evolve with it. Small particles remain 
coupled to the gas. Larger particles wit h sizes of 1 cm or larger (depending on the local 


properties of the gas) are uncoupled (e.g.. lAdachi et al.lll976l:IWeidenschilling||l977al:lRafikov 






















































23 


2004; 

Birnstiel et ah 

2 C 

10 

). 

(e.g., 

Windmark et ah 

2012 


Collisional processes enable some particles to grow to larger sizes 


2 OI 2 I: iGaraud et al.ll2013l) 


For large particles, interactions with the gas produce a radial drift and circularize the or¬ 
bits. The characteristic time sc ale for these interactions is f Adachi et alll976 : Weidenschilhng 
1977a : Chiang fc Youdin 2010 '): 


T. 


15 


1 km 


1 AU 


10 ® g cm ^ 


X. 


(51) 


where r is the radius of the particle, pg is the local gas volume density (normalized to the 
typical density of a minimum mass solar nebula), and 7^ is the orbital period of the local 
guiding center. For all but very large objects with r > 100 km, this time scale is small 
compared with the precession time. Thus, small particles uncoupled from the gas likely find 
most circular orbits. 


Numerical simulations of ensembles of small particles interacting with the gas sug¬ 
gest the surface density of the solids is somewhat steeper than the surface density of the 


gas (e.g.. 

Youdin & Chiane 

bO 

0 

0 

Brauer et ah 

2008; 

Birnstiel et ah 

2010 , 

2012 ; 

Laibe et al. 

2012 ; 

Pinte & Laibe 

2014) 

Here we adopt a standard power law: 


S(a) = Sof-) (52) 

\aoJ 

where Sq = 10 g/cm^. To derive basic time scales, we assume that particles in the disk have 
radii of r = 1 km and a density of 3 g/cm^. If they are stirred to their escape velocity. 


Vesc 


m/r = (87rGp/3) r 


130 


P 

1/2 

- ^ - 

3 g/cm^ 


.1 km. 


cm/s. 


(53) 


then we can use standard kinetic theory to estimate the collision time as {nav)~^ where n 
is the number density, a is the collisional cross-section, and v is the relative velocity. The 
damping time from collision is then 


T. 


damp 


prX^ 

2\/27rS 


3400 


(54) 


P 


- - 


■ a ■ 

6 

1 

M 

0 

-1 

■ M ' 

3g/cm^ 


-1km. 


UAuJ 


10 g/cm^ 


_lMo_ 


1 - 1/2 


yr. 


where we assume that collisions are inelastic and the disk scale height is proportional to 
the relative speed v, divided by the orbital freq uency. This expression ignores gravitation al 
focusing, which reduces the damping time fe.g.. IWetherill fc Steward Il993l: IOhtsukilll999l) . 
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While the damping time is only linearly dependent on particle size, the range of sizes can 
be considerable in an evolving planet-forming disk. For pebbles (r 1 cm), the damping time 
is very fast: Tjamp is formally shorter than a dynamical time. For large particles (>1 km), 
the time scale for collisional damping is fairly long, approaching the precession time scale 
for the gas when r > 10 km. However, these larger pa rticles interact with smaller particles 


throu gh dynamical friction and viscous stirring (e.g., lOhtsuki et al.N2002l: iGoldreich et ah 


20041 ). For these processes, the typical damping time for large particles is a factor of 10-100 
smaller than suggested by Equation fl54p . Thus, damping times for particles with r < 100 km 
are short compared with the precession time. 

For comparison, the gravity of the particle ring induces the binary to process with a 
period 


Tbir,- 


bin—pre 


4.0 X 10® 


r So 1 

-1 

^bin 

1/2 

Ct'lYl 

10 g/civ? 


Lo.2AuJ 


2.5 Q.bin_ 


yr [particle disk]. (55) 


We conclude that in a particle disk, as in a gas disk, conditions allow for damping to most 
circular orbits, so they may serve as the equivalent to circular Keplerian orbits. 

Eventually, the largest particles contain more than half the mass of all the solids in an 
annulus of the disk. At this point, damping t imes become much longer than stirring times; 


orbit a l evolution then becomes chaotic (e.g., [Chambers fc Wetherilll Il998l: IGoldreich et ah 


20041: iKenvon fc BromlevI 1200611 . Around a single star, chaotic systems eventually settle 
down into stable multi-planet systems. Aside from the impact of the inner unstable region 
and resonances, we expect a similar evolution of chaotic systems around binary stars (e.g., 
Quintana fc Lissauerl 120061) . 


4.3. Perturbations from massive planets 


In a particle disk, a circumbinary Jupiter-mass planet dramatically shortens the pre¬ 
cession time for the binary. If planetesimals orbit within a few times the binary separation 
and the massive planet orbits well outside this region (e.g., as in Figure [6l Obm ~ 0.2 AU, 
satellite at a ~ 1 AU and “Jupiter” at 2 AU), then the binary precession rate is 


d^bin ~ 


3 nij Obin^ 
4M a? 


(56) 


where Mj and oj are the planet’s mass and orbital distance. Substituting parameter values 
comparable to observed binary systems and consistent with the set-up in Figure EJ the 
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precession period (27r/'d7bin) is 
Tbin-pre ^ 1.2 X 10® 


■ M ■ 

1/2 

M 

-1 

®bin 

-1/2 

Oj 

_1M0_ 


_1 M Jupiter, 


[o.2AuJ 


L2AuJ 


yr. 


(57) 


Any binary precession from the gas giant pertnrber is a small pertnrbative effect that pre¬ 
serves the apsidal alignment between binary and planetesimal orbits. 


4.4. Planet formation in circumbinary disks 


Our analysis establishes time scales for gas and particles to damp to most circular orbits. 
With no precession of the binary, these orbits do not process. Although binary precession 
induces precession in most circular orbits, these orbits remained apsidally aligned with the 
binary. In this way, most circular orbi ts provide reference frames in which the dynamics of 
planet formation takes place (see also iPichardo et al.l 1200511200811 . For example, collisional 
damping and self-stirring of planetesima ls modihes the free eccentricity of circumbinary 


planetesimals (e.g., as in IWetherilll 1198011 : these processes have no impact on the forced 


eccentricity (and high-freq uency modes) driven b y the central binary. Gravitational stirring 


from distant planets fe.g.. IWeidenschillind Il989l) also excites eccent ricity, driving epicyclic 


motion in a manner similar to the binary itself flHeppenheimerl ll978L see Figure [6]). 


When ensembles of orbiting planetesimals attempt to grow into a planetary system, the 
relative velocity between particles, 

V ~ Cfree^K) (58) 

is a key parameter which establish e s the efficiency of gravitational focusing and collision out 


1 Stewart & Wetherill 

1988; 

Wetherill & Stewart 

1993; 

Ohtsuki 

1999; 

Ohtsuki et al. 


20021). Small relative velocities favor growth by mergers; large relative velocities favor de¬ 
struction. To estimate the boundary between these regimes, we rely on the specihc collision 
energy Q*^ required to disperse half the mass of a colliding pair of planetesimals to inhnity. 
For rocky material. 


10 ® 


1 -0.4 


1 km. 


+ 7 X 10® 


P 


- j. - 

3 g/cm^ 


-1 km. 


1 1.35 


erg/g 


(59) 


flDavis et al. 


1985 


Wetherill fc StewartI Il993l: iBenz fc Asphaugl Il999l: iHousen &: Holsapple 


19991: iKenvon &: Bromlevi 120051) . Equating to the specihc rest-frame collision energy 
between two equal-mass objects, mass loss exceeds mass gain when 


West ~ 0.1 km/s [destructive collisions, r = 1km] 


(60) 
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for kilometer-size rocky planetesimals. This relative speed is best measured in the reference 
frames of most circular orbits. 


Around a binary at orbital distance a, all planetesimals have an additional component 
to their velocity compared with their counterparts around single stars. This velocity, from 
the forced eccentricity, has a typical magnitude 


V "f^force — f'force^i, 


0.72 


{Mp- Ms) Ubin 


M 




1 


■ M ■ 

-1/2 

^bin 


6 bin 


■ a ' 

0.5 Mq 


_1 M©_ 


Lo.2AuJ 


Lo.2 J 


GAuJ 


km/s . 


(61) 

(62) 


Comparing this value with the disruption speed, it is clear that interpreting Uforce as a random 
motion leads to the simple prediction of destructive collisions for a broad range of particle 
sizes. Yet Uforce is the speed of a reference frame tied to a most circular orbit. Particles 
traveling on most circular orbits have no relative radial velocity; planetesimals may perturb 
each other and collide, but at velocities much smaller than Udest and Uforce, promoting mergers 
instead of destruction. 


This picture differs from the approach often taken in studies of circ u mbinary planet for- 


mation. For example, in 

Moriwaki & Nakaeawa 

(2004) 

. Meschiari 

12012) 

. Paardekooner et ah 

(2012 

) and 

Lines et ah 

(2014 

), particles in n-body simulations are initialized on circular Ke- 


plerian orbits about the binary center of mass. Gas, if present, is assumed to have fluid 
elements on exactly circular sub-Keplerian orbits. From ^ a particle trajectory with total 
eccentricity of e = 0 is identical to a trajectory with equal parts free and forced eccentricity 
(cfree = Cforce), where the phase of the free part is chosen to yield an initial net eccentricity 
of zero {'ipe = tt). The difficulty with these initial conditions is that the relative velocities 
are set with v ~ ~ 'yforce- It is then just a matter of time before precession of the free 

epicyclic motion drifts from the force motion, and particle orbits can cross. That time is 


T = 

^ pre 


27r 


67 


w 


flbin 

-2 

■ a ■ 

7/2 

■ M ■ 

1/2 


[o.2AuJ 


a AuJ 


_1Mq_ 


_0.5Mo_ 


-1 


yr, 


(63) 


where the reduced mass is /r < 0.5M. Entrainment of solids by the gas on snb-Keplerian or¬ 
bits c an reduce relative velocities, but only for particles of similar size fe.g.. lMarzari fc Scholl 


2000l ). Thus, setting up protoplanetary disks with all objects having comparable free and 


forced eccentricities quickly dooms them to destruction. 


Sometimes, the initial orbits of circumbinary planetesimals have no impact on the out¬ 
come of planet formation calculations. For example, if the free eccentricity is much larger 
than the forced eccentricity, then the binary’s time varying potential is an unnoticed pertur¬ 
bation. This situation applies in the late stages of planet formation, when most of the mass is 
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concentrated into a small number of large objects. At this stage, larg e objects have random 


velocities comparable to t he escape velocity of the largest object fe.g.. iGoldreich et ahll2004 


Kenvon &: BromlevI 120061) . Setting Uesc ~'i^force, we obtain 


Vo ^ 1000 


p 

-1/2 

■ M ■ 

1/2 

1 

1_ 


®bin 


6 bin 


■ a ' 

3 g/cm^ 


_lMo_ 


0.5 Mq 


Lo.2AuJ 


Lo.2 J 


CAuJ 


- 3/2 


km. (64) 


Outside of the critical radius for stability, Ocrit, a simulation with a collection of Pluto-size or 
larger objects should fail to see a difference between cir cumstellar and circumbinary environ¬ 
ments. The simulations of lQuintana fc Lissaueii (120061) support this interpretation. Except 
for the effects of the destabilization of orbits near Ocrit, ensembles of roughly lunar mass 
objects grow into a few terrestrial planets in circumstellar and circumbinary environments. 


4.5. Additional issues: resonances and binary evolution 


Once planetesimals have settled around most circular paths to grow by coagulation, 
the binary still i nfluences their evolution. Overlapping resonances may stir particles and 
destabilize them (1Wisdomlll980l) . However, these unstable orbits are typically close to the 
binary. For example, the 3:1 resonance, which corresponds to a singularity in the denomina¬ 
tor of Equation fl25|l . lies within 2.1a hin. As a guide, thes e unsta ble orbits are usually within 
the critical radius Ocrit identihed by iHolman &: WiegertI (1199911. In some binary configura¬ 


tions instabilities extend beyond Ocrit (e.g., iMusielak et al.l 120051: IPichardo et al.ll2005l. 12008 


Doolin fc Rlnndellll2011ll. Numerical simulations provide a straigh tforward way to identify 
these orbits (e.g. Popova fc Shevchenkol 2013 : Chavez et ah 201^ and to understa nd any 
impact on gas or solid particles in the planetary disk (e.g., IPierens fc Nelsonll2008bl) . 


Tidal evolution of the binary also complicates the long-term stability of circumbinary 
planetary systems. Changes in the orbital separation and eccentricity can change the po¬ 
sitions of critical resonances, leading to instability in systems which had been stable. In 
the Pluto-Charon system, satellites become unstable when tidal evolution drives a relatively 


rapid expansion of the binary orbit flWard fc Canupll2006l: Ihithwick fc Wull2008l: ICheng et ah 


20141) . For binaries composed of solar-type stars, tidal evolution on time scales of 1-100 Myr 
is im portant only for systems with orbital periods of 10 days or less flMeibom &: Mathieu 
20051 ). Thus, tidal evolution of the central binary probably has little impact on the formation 
of most circumbinary planetary systems. 
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4.6. Summary 


Unless the collapse of a molecular cloud into a binary + disk system is significantly 
different from the collapse to a single star + disk system, we expect the stages of circumbi- 
nary planet formation to parallel those of circumstellar planet formation. As long as non- 
axisymmetric structures such as spiral density waves are relatively unimportant, the time 
scales for gas damping, gas drag, collisional damping, and dynamical friction are all much 
shorter than typical precession times for the binary. Thus, material in a circumbinary disk 
should hnd its way onto most circular orbits around the binary in a similar way as material 
in a circumstellar disk damps onto circular orbits around a single star. Once material lies 
on or close to most circular orbits, small relative velocities of particles on adjacent orbits 
strongly favors growth over collisional disruption. 

Precession of the inner binary has little impact on this conclusion. When a massive 
disk or a giant planet causes the inner binary to process, the orbits of solid particles on most 
circular orbits maintain apsidal alignment. Because relevant damping time scales are short 
(< 10 yr for gas and < 10^ yr for particles) compared with the binary precession time scales 
(> 10^ yr for precession induced by a massive disk, or longer with particle disks or gas giants 
beyond 1 AU) both gas and solids can maintain this apsidal alignment (see Equations fITTI) 
and (IM)) for damping times and Equations flSUD . and flFTD for binary precession periods). 
Most circular orbits remain the reference frame in which to measure relative velocities. Thus, 
planet formation proceeds in a standard fashion. 


For calculations of circumbinary planet formation, starting with the right initial condi¬ 
tions for orbits of gaseous material or solid particles is crucial. Around a single star, it is 
sufficient to start material on orbits with modest eccentricity: precession is relatively unim¬ 
portant and all damping processes lead to circular orbits. Around a binary, it is important 
to differentiate between the forced eccentricity driven by the b inary and the free ec c entric - 
ity relative to a most circular orbit around the binary (see also iPichardo et al.ll2005l. l2008l) . 
Starting particles on circular Keplerian orbits around the binary center of mass creates a free 
eccentricity relative to a most circular orbit. In general, this free eccentricity is comparable 
to the forced eccentricity induced by the binary. For eccentric binaries, this extra eccentricity 
produces spuriously large collision velocities which can lead to collisional destruction instead 
of growth by merger. To provide a proper evaluation of circumbinary planet formation, it is 
essential to begin with most circular orbits and then evaluate relative collision velocities for 
appropriate values of the free eccentricity. 
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5. Application: The Kepler circumbinary planets 


In the last decade, data from the Kepler satellite have st arted to paint a rich picture of 


circumbinary pl anetary systems (e . g.. lArmstrong et a. 


known systems flDoyle et al.ll201ll: IWelsh et ah 


2 OI 3 I: iKostoy et al 


2013 


2012 


2014 and referen ce s therein). The 


Orosz et al.l l2012al]bl: ISchwamb et ah 


20141) have fairly massive planets with radii r ^ 0.25-0.76 Rj and 


semimajor axes a ~ 0.3-1.1 AU orbiting binaries with a broad range of mass ratios (0.25-1) 
and eccentricities (0.02-0.52). Although the number of circumbinary planetary systems is 
still rather small (~ 10), d etection rates suggest these planets are roughly as common as 
planets around single stars flArmstrong et al.l 120141) . 


At larger semimajor axes, circumbinary debris disks are also common. Among main 
sequence star s with FGK spectra l types, single stars and binaries are equally likely to haye 


often co-planar (e.g.. 


Kennedy et al. 

2012; 

Kennedy 

2015) 


Small number statistics cur¬ 


rently preyents robust conclusions, but it seems plausible th at this structure is primordial: 
circu mbinary disks form in the plane of the binary system (IKennedy et al.l l2012l: [Kennedy 

2OI5I ). 


Understandin g the structure of p rimordial circumbinary disks is also hampered by small 
sample sizes (e.g., iHarris et al.N2012h . Among multiple stars with ages of 1-3 Myr in the 
Taurus-Auriga molecular cloud, disks around single stars haye masses similar to those of very 
wide binaries wit h a > 300 AU (0.001-0.1 M c^, on the basis of dust mass estimates from 
millimeter fluxes; Andrews &: Williams 2005). Binaries with a ~ 30-300 AU (5-30 AU) 
haye factor of 3-5 (5-10) smaller disk masses. In three binaries with separations ranging 
from about 0.05 AU to 10 AU, disk mass estimates yield 0.01-0.03 Mq, comparable to the 
disks in single stars and wide binaries. For reference, a mass of 0.01 Mq corresponds to the 
Minimum-Mass Solar Nebula (MMS N), which is a the lower limit on the m ass required to 
build the planets in the solar system fjWeidenschillind Il977bl: iHayashil Il98ll) . 


Here, we suryey known Kepler circumbinary systems to interpret how their orbital 
characteristic s and size mi g ht inf orm us of their origin. While we reiterate some of the 
discussion in Leung fc Lee fj2013 ). we also focus on how the obseryations might impact 
our understanding of planet formation scenarios, whether these planets formed in situ or 
migrated from some larger semimajor axis inward. 


To frame the problem, we consider several different formation scenarios: 


1. In situ formation with no migration. Planets grow by coagulation from nearby gas and 
dust. The hnal planetary mass is limited by the initial surface density of the proto- 
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planetary disk. A Minimum-Mass Solar Nebula with Sq = 7 g/cm^ in Equation fimi 
has 4 Mffi in solids inside of the snow line a 2.7 AU and roughly 50 in solids inside 
50 AU fjWeidenschillind Il977bl : iHavashil llOSlIl . The mass in gas is roughly 800 
inside 2.7 AU and roughly 3000 inside 50 AU. Numerical simulations suggest rocky 
planet formation is inefficient; collisional processes often lose sig nificant amounts of 
mass during the assembly of Earth-ma ss and larger planets (e.g., iKenvon fc Bromlev 


2004 I2OO6I: [Raymond et al.ll2006l. I2OIII ). Augmenting the MMSN by a factor of 2-3 is 


sufficient to yield Earth-mass planets at 0.7-1 AU in the solar system. 

Producing a Neptune-mass planet in situ at 1 AU is more challenging. Accreting 
sufficient gas frorn the disk to produce Neptune requires a ~ 10 M m solid core (e.g., 


Pollack et al.lll996l: lRafikovll201ll: [Rogers et al.ll201ll: [Piso et al.[[2015[) . In the standard 


MMSN, the mass in solids is insufficient to produce such a massive core. When the 
surface density of the di sk is roughly 20 times the MMSN, coagulat ion models routinely 


yield 10 M 0 cores (e.g.. [Hansen fc Murravl [20121: [Schlichting[ [2014 see below). 


In most numerical simulations, formation of 2-3 Earth mass or larger planets is com¬ 


mon; producing a single planet is rare (e.g., [Kenvon fc Bromled [20061: [Raymond et ah 


201l[: [Hansen fc Murravl [2012[: [Hansen] [20141 1. If these simulations are ‘missing’ an im¬ 


portant piece of physics which allows several Earth-mass planets to merge into a single 
super-Earth mass planet, in situ formation of single super-Earth or Neptune mass plan¬ 
ets might be possible in lower mass disks, e.g., 3-5 times the mass of the MMSN instead 
of ~ 20 times. Without a better understanding, we assume that formation of Neptune 
mass planets requires a very massive disk. 


H. Migration then assembly. Precursor solid material is first moved fro m large a to within 
1 AU of the host star. Planet formation then proceeds in situ (iHansen fc Murray 


201211 . In this way, the limited amount of solids available inside of 1 AU is enhanced 


at early times; the delivery mechanism is uncertain. Because of the structure of the 
protoplanetary disk at small radii, gas accretion may be less efficient than beyond 
the snow line. Planets formed in this way may be “gas-starved” compared with more 
distant gas giants. 

HI. Migration through a gas disk. Gas giar its form beyond the snow l ine, where the solid-to¬ 
gas ratio is a factor of ~ 3 larger (e.g.- lKerinedv fc Kenvon[[2008]) . Planets move inward 
by exchanging torque with the gas disk flWardI [199711 . Composition and structure of 
these planets are typical of the solar system’s gas giants. 


IV. Planet-planet scattering. Massive planets formed beyond the snow line gravitationally 


scatter one another, either inward toward the central mass or outward ( 

Pasio & Ford 

1996; 

Chatteriee et ah 

2008; 

Juric & Tremaine 

2008 

Raymond et al.|[2010 

Marzari et ah 
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2010 ; 

Beauee & Nesvornv 

2012 ; 

Moeckel & Armitase 

2012 ; 

Bromlev & Kenvon 

2014) 


In a scattering event between two planets, if the more massive body is Neptune-size or 
larger (assuming typical solar system densities and orbital parameters for gas giants), 
then a broad range of outcomes is possible, including ejection from the system or the 
placement of a “hot Jupiter” near a single central star. High orbital eccentricity is the 
red flag for these scattering events. 


For each of these modes, the i nitial mass of the disk an d the epoch of planet formation 

and r eferences therein). I n the 


are important considerations (e.g., iNaiita fc Kenvon 


2014 


Taurus-Auriga molecular cloud (e.g., Andrews fc WilliamsI 120051: Andrews et al.ll2013l ). the 


median mass in solids for disks around single stars drops from 50-100 for protostars 
(ages of 0.1-0.5 Myr) to 10-20 M 0 for T Tauri stars (ages of 1-3 Myr). Observations of 


other star-forming regio ns suggest this evolution is typical (e.g., Andrews &: WilliamsI 120071: 


Williams &: Ciezal l201lh . Thus, existence of Neptune-mass to Jupiter mass planets favors 


early formation in the massive disks of young protostar^, 


5.1. Disk mass requirements for in situ growth 


For in situ formation, we estimate the surface density required to make a planet with 
mass m and radius r. Each protoplanet accretes material from a “feeding zone” with radial 
width 6a. Random motions of particles near the protoplanet set this width. When these 
particles have relative speeds smaller than the planet’s escape vel ocity, the planet c an accrete 
them. For particles with a free eccentricity, Cfree ^ ^^esc/w A.g.. ISchlichtind 120141) . Thus, 


Aa ~ 2aefree ~ 2a^ ^ 2a ( . (65) 

\ rM J 

where m and r are the core’s mass and radiu^. By choosing a form for the surface density 
(we use S ~ a~^'^ as in Equation 0371) 1 and integrating over the feeding zone (excluding any 
unstable region inside of Ocrit) we evaluate the minimum value for S at the planet’s position. 

Figure [8] shows the results of applying this prescription to the Kepler circumbinary 
planets in Table [U The minimum surface density required to build gas-accreting cores is 


"^For the solar system, radiometri c analyses of meteorites similarly suggests format ion of solids when the 
Sun had an age of 0.1-0.3 Myr fe.g.. iKleine et al.ll2009t iDauohas fc ChaussidoDll201ll) . 

®This result for the width of the feeding zone is so mewhat larger and probably more realistic than the 
width derived for the more standard “isolation mass” ( Greenzweig &: Lissauer lioonl) . 
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roughly an order of magnitude more than the MMSN. Thus, any in situ formation model 
requires a massive disk. Models which form these planets well outside the snow line require 
much less massive disks. 


This result allows us to eliminate in situ (no migration) models for all of the Keyler cir - 
cumbinary planets. If the Kepler circumbinary planets are common flArm.strong et al.ll2014f) . 
the precursor disks are also common. Among the youngests stars with ages of 0.1-0.5 Myr, 
however, disks with initial s urface densities of 10-20 times the MMSN are exceedingly rare 
flAndrews fc Williamsll2005h . Because lower mass protostellar disks are common, the Kepler 
planets require a formation model with some form of migration or scattering. 


5.2. Resonances, instabilities and migration 


Although formation followed by radial migration is a plausible path for the Kepler 
circumbinary planets, unstabl e orbital resonances in the circumbinary environment pose 
clear obstacles to migration fjPierens fc Nelson! 120071). Near certain c ommensurabilities , 
over lapping resonances e xcite large eccentricities fjWisdoiru Il982l. Il983l: iLecar et al.ll200lh . 
The iHolman fc WiegertI fjl999[) condition that orbits are unstable if they are within a dis¬ 
tance of Ocrit from the binary center of mass is a working guide; refinements are needed 
to identify the presence and “strength” of uns table resonances even outside of Ocrit (e.g., 
Popova fc Shevchenko 2013 : Chavez et ah 2015). 


Here, we assess whether the region around each binary in Table [T] is stable as a planet 
migrates inward. We consider two approaches: directly migrating the planet radially inward 
and, equivalently, expanding the binary. We show results for binary expansion models, which 
are more straightforward (and stringent) since they require adjusting only the Keplerian 
semimajor axis of the binary, instead of modifying the non-Keplerian orbit of the planet 
to mimic radial drift. Figure [9] shows the results for each Kepler planet at an equivalent 
radial drift rate of less than 1 0~^ AU/yr (typical for planets undergoing type I migration; 


Goldreich fc Tremainel Il980l: IWardI 119971: iTanaka et al.l 120021) The particles in the Figure 
manage to get as close to their binary host as the planet’s current position. Once inside, 
they become unstable near the 5:1 resonance, except for Kepler-34b, which becomes unstable 
near the 7:1 commensurability. Kepler-47b, orbiting the binary with the lowest eccentricity, 
is stable down to the 4:1 resonance. Thus, it seems plausible that the Kepler circumbinary 
planets, or smaller precursors (with fast damping times), could migrate through the potential 
mineheld of resonances. Further tests would be needed to conhrm whether the stability that 
we observe would remain with lower damping rates or slower migration times. 
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5.3. The Kepler circumbinary planets 

With constraints on in situ formation and migration from large distances established, 
we consider whether any of the fonr modes of planet formation is consistent with the bnlk 
properties and orbit of each Kepler circnmbinary planet listed in Table [T] We start with the 
hrst and most famons discovery, Kepler-16b. 


Kepler-16. The central binary, with a K-dwarf primary and an M-dwarf secondary, 
has significant eccentricity (cbin ~ 0.16). The planet is a gas giant with rad ins of 0.7 R, 
and a mass of 0.3 Mj, orbiting at 0.7 A U with low ec c entric ity, ~ 0.01 flDovle et ah 


201ll ). From the analytical theory (§2]), Ifvenng fc T^eel (l2013l) estimate that free and 


forced eccentricities are comparable, Cfree ~ ^toice ~ 0.03. 

The high mass and low eccentricity of the planet favor the migrate-then-assemble 
and the migrate-in-gas modes of planet formation. In situ formation with no migra¬ 
tion, althongh plansible, reqnires a very massive disk which is very nncommon among 
low mass stars with ages of 0.1-0.5 Myr (see Fignre |8]). Planet -planet scattering is 


less likely to resnlt in th is low-eccentricity confignration fe.g.. iMarzari et al.l 12010 


Beauge fc Nesvornvl 120121) . 


• Kepler-34. The binary has high eccentricity (cbm ~ 0.5) bnt the mass difference 
between the partners is small, within a few percent. Tims, the forced eccentricity at 
Kepler-34b’s location is low, aronnd 0.002. In contrast, the eccentricity of the planet is 
mnch larger, exceeding 0.2. The planet is abont the same size as Kepler-16b and may 
be a Satnrn analog. The high eccentricity and high mass favor a scattering scenario. 


• Kepler-35. This binary has stars that are also nearly eqnal in mass. The planet, 
Kepler-35b, is nearly the same size as Kepler-16b and Kepler-34b, bnt has a mnch 
smaller eccentricity (0.048). As with Kepler-16b, the migrate-then-assemble or migrate- 
in-gas modes are strongly favored over the scattering or in situ with no migration 
formation modes. 


• Kepler-38. The mass ratio of the binary is almost hve to one; the eccentricity is 
aronnd 0.1. The planet is comparatively small, with a radins of 0.4 Rj, snggesting a 
mass of ronghly 20 M®. Despite its low eccentricity (consistent with zero) and low 
mass, in situ formation with no migration is still problematic. Kepler-38b reqnires 
a snrface density that is more than a factor of twenty larger than the MMSN. The 
low eccentricity similarly eliminates the scattering mode. Thns, we favor either the 
migrate-then-assemble or migrate-in-gas mode. 
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Kepler-47. The binary consists of a Sun-like star and a red dwarf on a nearly circular 
orbit. The forced eccentricity of its planet is low; the eccentricity of the planet is 
consistent with zero. With a radius just under 3 R®, Kepler-47b probably has the 
smallest mass of all the known circumbinary planets, perhaps 10 M^. This planet 
has a companion about twice its radius, Kepler-47c, at an orbital distance of about 
1 AU. The combined mass of the two planets works against in situ formation with no 
migration. T wo planets orbiting inside the snow line requires a complicated scattering 


scenario see 


Moeckel fc Armitagell2012l) . The migrate-then-assemble and the migrate- 


in-gas modes are less complicated. Improved statistics for this kind of circumbinary 
planetary system would provide better constraints on the scattering mode. 


• PHl/Kepler-64. This binary has both a signihcant mass difference between its 
partners and a modestly high eccentricity (0.2). The resulting forced eccentricity is 
the largest of this sample (0.044). The planet PHlb is half the size of Jupiter, with an 
eccentricity of 0.05; thus, the free eccentricity is probably low. As with the other Kepler 
circumbinary planets, the large mass of the planet precludes in situ formation with no 
migration; the small orbital eccentricity eliminates most scattering models. Early 
formation in a massive disk allows either the migrate-then-assemble or the migrate-in- 
gas models. 


• Kepler-413. This low-eccentricity binary hosts a modest size planet, a third the 
radius of Jupiter. Like Kepler-38b, this mass is uncomfortably high for in szfu formation 
without migration. Unlike Kepler-38b, Kepler-413b has a signihcant free eccentricity of 
more than 0.1. Large eccentricity tends to preclude the migrate-in-gas mode. However, 
either the migrate-then-assemble mode or a scattering model can produce a planet 
similar to Kepler-413b. 


5.4. Summary 

The circumbinary planets observed by Kepler are Neptune-size or bigger and located 
just beyond the critical radius around their host stars. Although in situ planet formation 
with no migration is a promising way to grow Earth mass planets at these distances, it 
is very unlikely to be responsible for the known Kepler circumbinary planets (Figure [8]). 
The central issue — not enough mass to build Neptune-size planets — is circumvented by 
importing solids from beyond the snow line. Other formation mechanisms accomplish this 
mass transfer by invoking an inward radial hux of small particles within the disk, migrating 
fully-formed gas giants through the disk, or scattering gas giants from outside the snow-line 
where other large planets form. 
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Migration seems to be involved in most of the Kepler planets. However, withont larger 
samples of planets, it is impossible to distinguish between models where migration precedes 
assembly from those where migration follows assembly. All of the planets are too massive 
to allow in situ formation with no migration. However, the high free eccentricity observed 
in Kepler-34b and Kepler-413b are consistent with scattering events. Improved constraints 
on the orbits and bulk properties (mass, composition, spin, etc.) might allow more rigorous 
conclusions on their origin. 


_ ^Here our most important contribution to the discussion of close-in planets (e.g.,lSchhchting 

2 OI 4 J) is not to discriminate between formation mechanisms, but to emphasize that all are vi¬ 


able in the circumbinary environment at a > Ocrit- Thus, all issues concerning the formation 


of Neptunes and Jupiters inside the snow line for a single star carry over to the circumbinary 
case. 


6. Discussion 


Our main conclusion is that outside of a small region near a binary star, planet formation 
proceeds in much the same way as around a single star. This result stems from the existence 
of a family of “most circular” orbits around binaries that do not intersect, analogous to 
concentric circles around a single star. Gas, dust, and growing planets orbitally damp to 
these streamlined paths to avoid mutual collisions. The growth of planets, involving mergers, 
fragmentation, stirring and dynamical friction, takes place in the reference frame of guiding 
paths on these most circular orbits, just as it does in circumstellar disks. Without such 
paths, planetesimal orbits would inevitably mix at high velocity, leading to destruction, not 
growth. 


These rnost ci rcular paths are rooted in analytical theory. iLee fc Pealel (120061) and 


Leung fc Leel (120131) lay out the foundation, describing how circumbinary orbits are approxi¬ 


mated as linear combinations of (i) rapid, forced oscillations in response to the time-varying 
potential, (ii) slower, epicyclic motion — the “forced eccentricity” — that responds to the 
binary’s eccentricity, and (iii) “free eccentricity” and inclination relative to the plane of the 
binary. The periapse of the forced eccentric orbit is aligned with the periapse of the binary. 
The free eccentricity and inclination have the same meaning for a circumbinary orbit as 
in a Keplerian system. There is precession of the free eccentric orbit in the circumbinary 
case, for the same reason as in the case of a single, oblate star (quadrupole and higher order 
contributions to the time-averaged gravitational potential). 


We place this analytical framework in the context of planet formation. Our results 
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demonstrate that (i) most circular orbits exist for binaries even when their eccentricity is 
large enough that the theory — based on a perturbative approach — becomes questionable; 
(ii) non-intersecting, most circular orbits exist in the presence of a massive disk and a large 
planet; (hi) these orbits remain apsidally aligned with the binary — if the binary precesses, 
the forced epicyclic paths in the disk also process; and (iv) the response of a satellite to 
external perturbations (e.g., stirring from a distant planet) is the same as it would be around 
a single star, in the reference frame of a circular guiding center. Thus, planets can form in 
situ as close as a few times the orbital separation of their binary host s. Nearer to home, 
satell ites of the Pluto-Charon binary may have formed in a similar way flKenvon &: Bromlev 
2ni4h . 


Extending these ideas to include gas within a protoplanetary disk, disk gravity and gas 
pressure slightly distort the shape of most circular orbits but not their alignment with the 
binary. Alth ough deviations from strictly circular geometry introduce new hydrodynamical 
effects fe.g.. lOeilvie fc Barken I2014J: [Barker fc Ogilvid 1201411 . circumbinary gas streamlines 
plausibly settle on these most circular paths as they do on circular orbits around a point 
mass. Interactions with solid particles are also similar. The “headwind” felt by planetesimals 
traveling through the gas is nearly the same as in a circumstellar disk (Figure E]). As the 
gas dissipates in time, streamlines and planetesimals damp to the same set of most circular 
orbits. 

To explore whether our resul ts are sensitive to instabilities of circumbinary orbits, w e 


also consider resonant excitations flArtvmowicz fc Lubowlll994J: iKlev &: Haghighipouiil2014j) . 


In our tests, the Kepler circumbinary planets are all beyond the outermost unstable reso¬ 
nance around their host. For most systems, th is resonance is the 5:1 commensurability, close 


to the prediction of lHolman fc WiegertI fjl999l) . We dehne stability over a limited time frame 


with typical damping rates 


fWeidenschilling 

1977b: 

Goldreich et ah 

2004; 

Chiang & Youdin 


2010 ) and migration times ( Ward 1997I) . both in situ formation and migration of planets or 


their precursors are plausible. 

Issues that we do not address include possible non -axisymmetric structure in a massive 
circumbinary disk (jPelupessv fc Portegies Zwart 2013 ). or the physics near the edge of the 
stability zone at a distance of a few times the binary separation. Simulations of gas and 
particle disks can shed light on whether material is pushed outward by torque exchange wit h 


the binary or lost in the unstable zone f Artvmowicz fc Lubow 19961: Giinther fc Kiev 


Furthermore, we do not consider tidal evolution of the binary (e.g., iHurlev et ah 


2002 ). 


20021 ) . 


Depending on the nature of the evolution, variations in binary separation and eccentricity 
will affect the stability of orbits in circumbinary planetary system. For an intriguing example. 
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see the discussion bv lWard fc Canupl (120061) on resonant transport of moo ns in orbits around 
the Pluto-Charon binary as it tidally expands flBromlev &: Kenvonll2015l) . 


Our results stand in contrast to previous theoretical studies of circumbinary planet for- 


to destructive collisions 

(Moriwaki & Na 

caerawa 

2004; 

Meschiari 

2012 : 

Paardekooner et ah 

2 Q12; 

Rafikovibnbi 

Xie 

2Q13; 

Lines et ah 

2014 

). In this interpretation, if the binary’s eccen- 


tricity is even modestly high, Cbin ~ 0.2, epicyclic velocities eventually exceed the shattering 
speed for all but the largest planetesimals over a s ignihcant range of orbital distan ces. En- 
train r nent in a gas disk (e.g..lMarzari fc Schollll2000[) , or precession from disk gravity flRahkov 
20131: ISilsbee fc Rahkovi l2015a|) mitigate the situation only for some range of planetesimal 
sizes or distant regions of a protoplanetary disk. In this view, planet formation close to the 
binary is impossible. 

These conclusions depend on the assumption that gas and pl anetesimals are in itially on 
orbits with no eccentricity in the Keplerian sense (see Figure 1 in lLines et al.ll2014i) . In light 
of analytical theory (§2), this choice endows particles with equal amounts of free and forced 
eccentricity. As the free part drifts in phase, collisions destroy growing protoplanets. A more 
realistic approach to modeling an unstirred disk is to set the gas and planetesimals on most 
circular orbits with no free eccentricity. Absent any stirring, planetesimals in a particle disk 
orbit an eccentric binary host indehnitely without colliding. If planetesimals develop small e 
and i about the most circular orbits, collision velocities are modest, just as for planetesimals 
with modest e and i about a circular orbit around a single star. Modest collision velocities 
promote growth instead of destruction. Planet formation close to the binary is then robust. 

Finally, we review characteristics of Kepler circumbinary planets in light of our results. 
If circumbinary planets form roughly in situ, the orbital characteristics and sizes of these 
planets require very large initial surface densities (10-20 times the MMSN). Although these 
high surface densities ap pear to preclude in situ m odels with no migration, the migrate- 
then-assemble picture of iHansen &: Murray! (120121) is viable. As a plausible alternative, 
migrate-in-gas models allow these planets to form beyond the snow line and migrate inward 
through the gas. For this set of Kepler planets, migration through gas avoids unstable 
resonances around the binary. For Kepler-34b and Kepler-413b, formation beyond the snow 
line followed by a scattering event is also a reasonable scenario. 

In general, beyond the inner unstable cavity around a stellar binary, the evolution of 
solids on most circular orbits differs little from the evolution of solids on circular orbits around 
a single star. Thus, the standard issues of formation around single stars (e.g., planetesimal 
formation, migration, resonances, scattering, et c) have clear para llels in circumbinary disks 
(for a summary of the issues in single stars, see ISch]ichtingil2014j) . 
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We conclude with predictions for circumbinary planetary systems. In our scenario, 
planets are as prevalent around binaries as around single stars. Furthermore, relative to 
coplanar, most circular orbits, these planets should have the same distribution of orbital 
elements (free eccentricity and inclination) as their circumstellar cousins. Data from the full 
Kepler catal og indeed suggest that the planets have comparable rates around binaries and 
single stars flArmstrong et al.ll2014l) . Tatooine sunsets may be common after all. 
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Table 1: Kepler binaries. Along with orbital parameters are estimates of the innermost stable 
orbit radius (ocrit) and the forced eccentricity (cforce)- 



Mp (Mq) 

Ms (Mq) 

Q^bin (-^U) 

Cbin 

a(AU) 

e 

r(Rj) 

Q^crit (-^U) 

Cforce 

Kepler-16“ 

0.687 

0.202 

0.224 

0.160 

0.720 

0.024 

0.75 

0.646 

0.034 

Kepler-34“ 

1.049 

1.022 

0.228 

0.521 

1.086 

0.209 

0.76 

0.833 

0.002 

Kepler-35“ 

0.885 

0.808 

0.176 

0.142 

0.605 

0.048 

0.73 

0.496 

0.002 

Kepler-38'^ 

0.949 

0.249 

0.147 

0.103 

0.464 

<0.032 

0.39 

0.389 

0.024 

Kepler-47“ 

1.043 

0.362 

0.084 

0.023 

0.296 

<0.035 

0.27 

0.203 

0.004 

PHl'^ 

1.528 

0.378 

0.174 

0.212 

0.634 

0.054 

0.55 

0.527 

0.044 

Kepler-413® 

0.820 

0.542 

0.099 

0.037 

0.355 

0.118 

0.39 

0.253 

0.003 


“See Dovle et al. ( 20111 ) and I Welsh et all (j2012ri : orbital elements are from iLeung fc Leel (1201,11 Table 1); 
^ Orosz et al. (2012b). 


Orosz et al.l (l2012al) : A second planet (r = 0.41 Rj) is at ^ 1 AU. 


‘ Schwamb et al. ( 2013ll and Rostov et al.l ( 2013 ): Kepler-64. 


Rostov et all (120141) . 
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time (orbital period) 

Fig. 1.— The radial excursion of Kepler-16b on most circular orbits. The dark curve 
shows the radial excursion of a satellite at the orbital position of the planet in the absence 
of free eccentricity and inclination, plotted as a fnnction of orbital phase (in nnits of the 
planet’s orbital period). The gray cnrves show orbits at slightly displaced orbital distances. 
The trajectories are dominated by the forced eccentric orbit. Comparatively small higher- 
freqnency oscillations are visible in the curves. Despite the appearance that these oscillations 
are not in phase between the curves, the orbits do not cross. 
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Fig. 2.— Illustration of radial excursions of a satellite on a most circular orbit. The orbital 
configuration and binary masses are derived from Kepler-16 and its planet Kepler-16b at 
a = 0.70 AU. The binary eccentricity has been reduced by a factor of 10 from the real 
system; the plot then distinguishes the forced eccentricity (larger amplitude, driven at the 
orbital period, ~ 0.6 yr) and the high-frequency oscillations (smaller amplitude driven at 
the binary’s orbital period and the synodic period of the satellite relative to the binary, 
~ 0.1 yr). The black curve is from simulation while the gray curve is from analytical theory 
(Equation ([21])). 
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Fig. 3.— Comparison between a most circular orbit and one initially on the circular path of 
the guiding center (zero eccentricity). The most circular path (upper curve) is non-precessing, 
with excursions from the guiding center that have a maximum amplitude that does not drift 
over time. In contrast, a particle that is set up on a “circular” orbit — launching the particle 
from its guiding center with the speed of uniform circular motion about the binary center of 
mass — has a mixture of free and force eccentricities in equal measure. The relative phase 
allows zero eccentricity at the start. The result (lower curve, displaced from the upper one 
for clarity) is the beat pattern with a frequency given by the precession rate of the free 
eccentric orbit (Equation ([33])). 
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Fig. 4.— Simulated most circular orbits around a binary with moderately high eccentricity. 
The orbital conhguration is derived from Kepler-16, as in the previous hgures, but the binary 
eccentricity is set to 0.5, and the satellite is placed at 0.9 AU, just beyond the critical radius 
for stability (ocrit = 0.86 AU). The center of mass of the system is at the origin in this 
X — y map of the plane of the binary. The secondary’s pericenter is on the positive a:-axis. 
The satellite (blue dots) tracks a narrow elliptical path with an eccentricity of Cforce ~ 0.08, 
hxed and aligned with the binary for ~ 10^ orbital periods (~ 10^ binary orbits). When 
the potential of a massive disk (S of 2000 g/cm^ at 1 AU) is included, the satellite’s most 
circular path does not process (magenta points, mixed in with the blue ones). For reference, 
we show samples of a satellite orbiting a single star with Kepler-16’s total mass, also in 
this disk potential (gray points). In this case, the precession of the satellite’s argument of 
periastron is rapid (~ 0.005 rad/yr); the points £11 an annular swath. In these simulations, 
there is no interaction between the disk and the binary. 
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Fig. 5.— Inspiral from gas drag. A sub-Keplerian gas disk has fluid elements on most 
circular orbits with r] = 0.001 (Equation fHOli i. Embedded in it is a particle initially on a 
most circular orbit at 1 AU (outer blue ring). The particle evolves, drifting inward as a result 
of a drag force proportional to its speed relative to the gas (light blue points; the local gas 
speed is calculated using the Lee-Peale-Leung analytical theory). The particle’s hnal orbit 
(inner blue ring) remains apsidally aligned with the binary, a pair of stars like Kepler-16 
except with Cbm = 0.5 (inner black curves). A circular path (magenta curve) provides a 
reference. 
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Fig. 6.— The effect of an external pertnrber on circnmbinary orbits. The cnrves show 
most circnlar paths aronnd Kepler-16 at the orbital distance of Kepler-16b, both in isolation 
(magenta cnrve) and with a Jnpiter-mass planet orbiting with e ~ 0 at 2 AU (bine cnrve). 
The difference between these two trajectories (cyan curve; offset for comparison) compares 
well with data from a satellite orbiting a single star with the mass of Kepler-16 and a 
Jupiter-mass companion (orange curve). These lower curves show that the most circular path 
provides a frame of reference for the action of external perturbations, just like a circularly 
orbiting guiding center in the circumstellar case. 
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Fig. 7.— Binary precession and circumbinary orbits in a simulation. The outer curves show 
most circular paths (a = 0.9 AU) around a binary like Kepler-16, set with eccentricity of 0.5 
(inner curves) and with a Jupiter-mass planet orbiting with e ~ 0 at 2 AU. The gray and 
light blue-shaded curves show orbits in the inertial reference frame of the system’s center of 
mass. The duration of the simulation is about a quarter of the precession period; orbits have 
precessed about 90° (see Equation (l56|) h The black and dark blue ellipses are the same orbits 
represented in a reference frame that precesses with the binary’s periapse. The satellite’s 
forced epicyclic motion evidently precesses at this rate. 
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Fig. 8.— The minimum surface density to build the Kepler circumbinary planets. Each 
planet is shown at its orbital distance from the host binary. The value of the surface density 
(S) comes from determining an annular width from which the planet could have accreted 
mass, based on its escape velocity when it had the mass of a 10 M^, prior to its acquisition 
of a gas atmosphere (from Equation (l65|l ). The solid line is the surface density of a Minimum 
Mass Solar Nebula (S = 7(a/l g/cm^), while the dashed and dotted lines correspond 

to di sks with three and ten ti mes that density (as labeled). Indications from simulations 
fe.g.. lKenvon fc Bromlevll2006h suggest that the intermediate-mass disk is a realistic starting 


condition for the Solar nebula. 













